Singe ferroelectric and chiral magnetic domain of 
single-crystalline BiFeOs in an electric field 

Seongsu Lee/ Taekjib Choi/ W. Ratcliff II/ R. Erwin/ S-W. Cheong/ and V. Kiryukhin^ 

^Rutgers Center for Emergent Materials and Department of Physics and Astronomy, 
Rutgers University, Piscataway, NJ 08854, USA and 
"^NIST Center for Neutron Research, 
NIST, Gaithershurg, MD 20899, USA 
(Dated: September 11, 2008) 

Abstract 

We report polarized neutron scattering and piezoresponse force microscopy studies of millimeter- 
sized single crystals of multiferroic BiFeOa. The crystals, grown below the Curie temperature, 
consist of a single ferroelectric domain. Two unique electric polarization directions, as well as the 
populations of equivalent spiral magnetic domains, can be switched reversibly by an electric field. 
A ferroelectric monodomain with a single-g single-helicity spin spiral can be obtained. This level 
of control, so far unachievable in thin films, makes single-crystal BiFeOs a promising object for 
multiferroics research. 

PACS numbers: 77.80.-e,77.80.Dj,77.80.Fm,75.80.+q 
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Understanding and control of domain structure is crucial for applications of magnetic 
and ferroelectric materials, as well as for studies of their intrinsic properties. Coercivity of 
permanent magnets, switching of memory devices, and many other functional properties are 
defined by domain physics. Ferroelectricity and magnetism coexist in multiferroics[l]. These 
materials have attracted significant attention recently because of an intriguing possibility 
of control of their magnetic properties with an electric field and vice versa^ 2, y, 14] leading 
to discovery of new physical phenomena and fabrication of novel devices'S]. So far, little 
research has been done on the domain structure of multiferroics. Among a few notable 



exceptions 



6|] is BiFeOs (BFO), but mostly in its thin film form, since suitable single crystals 



were unavailable. BFO is arguably the most extensively studied multiferroic due to its large 
polarization at room temperaturej4]. Importantly, this is the only multiferroic exhibiting 
large effects of an electric field on the magnetic structure, as first demonstrated in thin 

poor crystallinity and are influenced by substrate strain, [2|, |5|] impeding studies of BFO's 
intrinsic properties, and complicating device fabrication. Herein, we report studies of single- 
crystal BFO in which an unprecedented level of control of domains and their electric and 
magnetic properties can be achieved. 

BFO is ferroelectric for T<Tc^850 °C. It exhibits the rhombohedral R3c structure with 
a perovskite pseudocubic unit cell (a~3.96 A, a^89.4°) elongated along the (111) direc- 
tion which coincides with the electric polarization vector P. Qj For T<T/v~370 °C, Fe spins 
exhibit the G-type antiferromagnetic order 10|. In single crystals, there is an additional 
long-range (A^620 A) magnetic modulation with the wavevector ti={6,0-6), believed to 
be cycloidal with antiparallel spins slowly rotating in the plane containing P and r.j^ 
Wavevectors T2={S-6,0) and T3={0-6,6) are equivalent by symmetry, giving three propaga- 
tion directions, each perpendicular to P. The spins are rigidly coupled to the crystal lattice, 
and rotation of P between the different cubic body diagonals by an electric field results in 
the rotation of the magnetic structure. ^, Q, Q] This effect was recently used to demonstrate 
local control of magnetization by an electric field, 

BFO exhibits four structural variants (ferroelastic domains) based on the four different 
cubic body diagonals, and eight possible directions of P (ferroelectric domains) shown in 
Fig. 1(a). To produce any realistic device, and to study the intrinsic properties, single- 
domain samples are needed. However, despite significant efforts, the best available BFO 
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thin films always consist of /im-sized ferroelastic domains. 0, 5| Films with suitable electrical 
characteristics are of the (001) type, and therefore poling (electric field applied along the 
cubic direction) is inappropriate for getting rid of the ferroelastic domains. In single crystals, 
poling in the (111) direction might be possible, but the situation is complicated by the 
presence of the equivalent magnetic domains in a single ferroelectric domain, as shown in 
Fig. 1(b). In this work, we report preparation of mm-sized single crystals with a single 
ferroelectric and single-r magnetic domain, see Fig. 1(c). The magnetic state is chiral, as in 
a single-helicity spiral. The electric and magnetic state can be reversibly controlled by an 
electric field, which switches the electric polarization between the two unique states shown 
in Fig 1(d), and also controls the equivalent magnetic domains. These results open new 
opportunities for studies of BFO, and for BFO-based device fabrication. 

Single crystals of BFO were grown using a Bi203/Fe203/B203 flux by slow cooling from 
870 to 620 °C. Crystals with ~mm^-area natural faces normal to the cubic directions were 
obtained. Piezoresponse force microscopy (PFM) was performed with a commercial scan- 
ning probe microscope. Px, Py, and Pz components of P were measured, with Pz normal 
to the sample surface, and the axes coinciding with the pseudocubic directions. For electric 
domain switching, a 40 V bias (bottom electrode) was applied to a 20 /im-thick sample, 
with the PFM tip grounded. Neutron experiments were done at T=50 K on BT9 and BT7 
(polarized measurements) triple-axis spectrometers at NIST Center for Neutron Research 
at an energy of 14.7 meV with collimations 40-10-S-40-80 (BT9) and 25-10-S-50 (BT7). The 
collimations were set such that their further tightening would result either in no significant 
resolution improvement, or in unacceptable signal losses, ensuring the best possible resolu- 
tion compatible with our signal level. Two parallel gold contacts were evaporated on the 
largest faces of a 3x2x1 mm^ sample for an electric field E to be applied in the (001) 
direction. 

Different ferroelastic domains can be easily detected using our neutron scattering setup 
[3]. The only Bragg peaks found in our samples were due to a single domain, even though 
a careful search for other domains was done. Considering the experimental accuracy, we 
conclude that more than 95% of each studied sample consist of a single ferroelastic domain. 
Multiple PFM measurements show the same P^, Py, and Pz across the surface in the majority 
of as-prepared samples. Typical images are shown in Fig. 2(a-c). Thus, the samples consist 
of a single ferroelectric domain. Local electric polarization can be switched reversibly by 
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application of E normal to the (001) surface by the PFM tip. Fig. 2(d-f) show that while 
Pz changes its direction in the region of bias application, and Py remain unchanged in 
the process. Thus, P rotates by 71° as shown in Figs. 1(d) and 2(g), the rotation direction 
is unique, and no 109° or 180° domains are created. This observation is in agreement with 
theoretical calculations showing that the elastic energy losses are minimal for the 71° domain 
boundaries. [11 1 



We now turn to the magnetic structure of the equivalent magnetic domains ri, T2, and 
r3 in the ferroelectric domain Pa of Fig. 1(b). Fig. 3 shows unpolarized neutron scattering 
pattern near the (1/2,-1/2,1/2) ±t magnetic peak in the {h,-h,l) plane, with projections of 
all the three magnetic wavevectors marked. Only the PA-domain signal is detected in this 
geometry, even when the second ferroelectric domain (Pb in Fig. 1) starts appearing for 
E<0. (See Ref. 7| for the discussion of the Pe-domain magnetism.) For E=0, two magnetic 
domains, with wavevectors ti and T2, are present. Application of E at an acute angle of 
55° to P, as shown in Fig. 1(c), suppresses domain T2. For E=1.3 MV/m, only domain 
Ti is observed within the accuracy of our measurement (~10% of the domain population), 
see Fig. 3(b). This figure shows that a fit to two Gaussian peaks at ±ri, with the widths 
and orientation closely matching the calculated resolution, reproduces the data very well. 
Note that for the out-of-scattering-plane vector ti, the peak positions differ slightly from 
their projections on the scattering plane due to resolution effects. Thus, the entire crystal 
consists of a single-r magnetic domain, as shown schematically in Fig. 1(c). Sufficiently 
large negative E produces the opposite effect, enriching domain T2 at the expense of ti, see 
Fig. 3(c). This effect is probably driven by the uniaxial piezoelectric strain which is linearly 
coupled to E. Figs. l(e,f) show that the field projection on the plain containing ti, T2, and 
T3 is normal to T2, and therefore the induced strain lifts the symmetry making the domains 
equivalent in zero field. Note that T2 is the only common magnetic wavevector for Pa and 
Pb domains, which may play a role in its stabilization for E<0. The absence of domain 
T3 is probably related to the post-growth residual strain. To characterize the magnetic 
domain populations, the neutron data were modeled by four Gaussian peaks of the same 
shape (measured at E=1.3 MV/m), see Fig. 3(a-c). For ri, the peaks were located in the 
experimental positions found at E=1.3 MV/m. For the in-plane T2, the peaks were located 
in the theoretical positions. As functions of E, the domain populations form reproducible 
hysteresis loops, showing electric field control of the magnetic domains, see Fig. 3(d). While 
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higher-resolution measurements are needed to reveal the details of these loops, our data 
strongly suggest that any specific magnetic domain can be stabilized by application of the 
appropriate uniaxial stress, and that complete control of the magnetic state can thus be 
achieved in BFO. 

I — I 

Cycloidal magnetic state is believed to be realized in BF0,|3] but alternatives were also 



proposed. 



12| While recent single- crystal neutron studies support the cycloid, [Sj non-coUinear 



states are difficult to analyze, and further work is needed to establish the magnetic state un- 
ambiguously. Various magnetic interactions can produce the cycloid. Common cases include 
competing nearest and next nearest neighbor exchange, and antisymmetric Dzyaloshinsky- 
Moria (DM) and magnetoelectric interactions. [1] In the (1-10) FeO chains in ferroelectric 
BFO, all the exchange- mediating oxygens are shifted from symmetric positions with respect 
to Fe in the same direction along P.jsl] For such shifts, the magnetoelectric interaction is 
proportional to P{myLx — rrixLy), where m is the unit magnetization, L is^ the antiferro- 
magnetic vector, and the xy-p\a.ne is defined as normal to the direction of P. [13] It is of the 
same form as the DM interaction, and should produce a weakly ferromagnetic spiral jl3. 
with rotation direction (helicity) uniquely defined in the rotation plane by P ~ e x r, 
where e is the spin rotation axis,[l| see insets in Fig. 4. Non-zero P is essential for the 
chiral state, i.e. the spiral is induced by P which appears at T much higher than T^. 
(Note that this is different from magnetically-driven ferroelectrics such as RMnOs, in which 
the situation is reversed, and P is induced by the inverse DM interaction in a magnetic 



spiral. 



ij) We used polarized neutron scattering to characterize the helicity of the spiral. 



Scans were taken in both of the spin-fiip channels, denoted as (+-) and (-+), with neutrons 
polarized parallel to the scattering vector Q. Scattering cross section[15[] is proportional to 
{S±YSly) + {S^zSlz)±^i{S±YSlz) - {S±zSly)) where {AB*) = E exp(zQrj)A(0)S*(rj), 
S± = Q X (S X Q), S is the spin at the position rj, and the coordinate axes X and Z run 
along Q and normal to the scattering plane, respectively. It depends on the helicity of the 
spiral and on the angle between e and Q: the difference between the channels is maximal 
for e||Q, and is absent for e ± Q. Fig. 4(a) shows scans through the (1/2,-1/2,1/2) ±ri 
peaks in the {h-h,l) plane in the single-domain state of Fig. 3(b). The direction of this scan 
is shown in Fig. 3(b) with dashed line. In this geometry, the intensity ratio of the +ti and 
-Ti peaks, calculated for the spiral rotating in the plane of P and ri, is ~8. The data of 
Fig. 4(a) fit well to the sum of these peaks with the intensity ratio of 9, and a weak central 
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peak, possibly coming from A/2 contamination. Note that for any non-chiral structure, such 
as a spin density wave, or for a sample equally populated with spirals of the both helicities, 
both the Ti peaks should be of the same intensity. Fig. 4(b) shows scans through a different 
magnetic peak, (1/2,1/2,1/2) ±r, in the {h,2l-h,l) scattering plane in the two-domain state 
of Fig. 3(a). As discussed in Ref. j^, the ti and T2 domains produce peaks in the ±6/2 
and ±6 positions, respectively, in these scans, giving rise to the broad peaks of Fig. 4(b). 
In this geometry, e is normal to Q, and no difference between the (+-) and (-+) channels is 
observed, as expected. Thus, our data provide strong evidence for the spiral model discussed 
in the literature, jo] Moreover, we show that magnetic domains are fully chiral (contain a 
single-helicity spiral for a given r), as expected on symmetry grounds. 

Spontaneous formation of large single crystals containing a single ferroelastic, ferroelec- 
tric, and (upon application of an electric field) single-r chiral magnetic domain, while seem- 
ingly surprising, was achieved by design. B2O3 was added to the flux to reduce the liquidus 
temperature down to ~620 °C, thus allowing the crystals to grow below Tc in the polar 
state, jl^ Effectively, the system behaves as a pyroelectric, and formation of a unipolar fer- 
roelectric domain is favored. While the growth temperature was not measured directly, we 
found that annealing the crystals at T^Tq produces ferroelastic domains, strongly support- 
ing this mechanism. The origin of the chiral magnetic state lies in the symmetry of the 
system, which fixes the helicity of the spiral once P and r are chosen, as observed in our 
experiments. 

The above data demonstrate that high level of control of electric and magnetic properties 
can be achieved in single- crystal BFO. The samples contain a single ferroelectric domain, 
reversible ferroelectric switching occurs between only two well defined polarization states, 
and a single-r magnetic chiral domain can be stabilized by piezoelectric strain. Note that 
50-/im-thick single-domain crystals were reported in Refs. |8|, [l6| , but no reversible control 
of P and no control of the equivalent magnetic domains was achieved, and the chirality was 
not measured. Also, as shown previously, the magnetic structure rotates together with P 

tin 

when an electric field is applied, p, |8| Such level of control is currently unachievable in thin 
films which consist of a patchwork of ferroelastic microdomains, with domain-dependent 
properties in an applied field. Single- crystal BFO is therefore suitable for studies of 
intrinsic properties of this multiferroic, as well as for fabrication of devices with precisely 
controlled properties. Our results may also help identify prospective approaches for growth 
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of uniform thin films. Finally, BFO provides an intriguing example of successful "chiral 



resolution" 



171], in which a large single crystal is prepared in a chiral magnetic state. This 



result reveals the symmetries involved in the coupling between the magnetic and electric 
order parameters in BFO. Chirality and the inversion symmetry are strongly coupled in 
BFO, as they are in many other multiferroics, such as TbMnOa, in which the chirality is 
also coupled to an electric field. 18(] Functional properties of these materials are defined 
by the underlying symmetry, and understanding the symmetry is essential for control and 
application of these properties. 

In summary, we report preparation of mm-sized single crystals of BFO which grow below 
the Curie temperature and consist of a single ferroelectric domain. The electric polarization 
can be switched by an electric field applied along the cubic direction, and only 71° rotation of 
P is observed. Magnetic state in BFO is chiral, and for a given wavevector the spiral is of a 
single helicity defined by the symmetry of the system via the direction of P. The equivalent 
spiral domains can be reversibly controlled by an electric field, and a chiral magnetic state 
with a single-r spiral can be stabilized by field-induced strain in the entire crystal. This is 
a highly unusual way of achieving chiral resolution and of controlling the magnetic state. 
These properties of single-crystal BFO make it a promising system for multiferroics research 
and experimental device fabrication. 

This work was supported by the NSF under Grants No. DMR-0704487, and DMR- 
0520471. 
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FIG. 1: (Color online) (a) Possible electric polarization directions in the pseudocubic cell of BiFeOs. 
(b) Equivalent magnetic spirals in a single ferroelectric domain, (c) A single ferroelectric/magnetic 
domain state, (d) Switching of the ferroelectric domains in an applied field, (e) and (f) show the 
orientation of E with respect to ri, r2, and rs, and its projection Ein on the plane containing 
these vectors. For simplicity, only one antiferromagnetic sublattice in the spirals is shown, and the 
modulation wavelength is decreased. 
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FIG. 2: (Color online) Representative AFM (a) and PFM (b,c) images of the (001) surface. Out- 
of-plane (d), and in-plane (e,f) PFM images of a ferroelectric domain created by a biased PFM 
tip. Arrows indicate the polarization direction probed in each image. The polarization vectors in 
the ferroelectric domains are shown schematically in (g). 



FIG. 3: (Golor online) Magnetic peaks near the (1/2,-1/2,1/2) position in the {h,-h,l) plane for 
E=0 in the as-prepared sample (a), for £^=1.3 MV/m (b), and for E=-1.3 MV/m (c). The ellipse 
in (b) shows the calculated instrumental resolution. Crosses indicate projections of the ri, T2, and 
Ts peaks on the scattering plane. Circles show the positions of the peaks used in the fits described 
in the text. The fits and the residuals are shown to the right of the main panels. Populations of 
the magnetic domains n and T2 determined in these fits are shown in (d) as functions of E. 



FIG. 4: (Color online) Polarized neutron scans in the (-1—) and (—1-) spin flip channels, (a) {h,-h,1.5- 
2h) scan through domain ri near the (1/2,-1/2,1/2) position in the {h,-h,l) plane. The direction 
of this scan is shown in Fig. 3(b) with dashed line. Error bars are from counting statistics, 
(b) (0.5-|-A/i,0.5-A/t,0.5) scan through domains ri and T2 near the (1/2,1/2,1/2) position in the 
{h,2l-h,l) plane. Arrows show the positions of the peaks due to these domains. Insets show the 
orientation of the scattering vector Q with respect to the plane of the magnetic spiral, as well as 
vectors e, P, and r. 
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